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Abstract 

Quantum phase transitions in a system of bosons with angular momentum L = 0, 2 
(s,d) and a single fermion with angular momentum j are investigated both classically and 
quantum mechanically. It is shown that the presence of the odd fermion strongly influences 
the location and nature of the phase transition, especially the critical value of the control 
parameter at which the phase transition occurs. Experimental evidence for the U(5)-SU(3) 
(spherical to axially-deformed) transition in odd-even nuclei is presented. 
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1. Introduction 

Quantum phase transitions (QPT) are qualitative changes in the structure of a physical 
system induced by a change in one or more parameters that appear in the quantum Hamilto- 
nian describing the system. Originally introduced in nuclear physics [H, [ij, where they were 
called ground state phase transitions, they have received recently considerable attention in 
condensed matter physics and other areas Quantum phase transitions were investigated 
both classically 0, [El and quantum mechanically in even-even nuclei in the early 1980 's 
within the framework of the Interacting Boson Model (IBM), a model of nuclei in terms 
of correlated pairs of nucleons with L = 0,2 treated as bosons {s,d bosons) f?\. In recent 
years, this study has been greatly expanded, including the aspect of symmetries (critical 
point symmetries fsl. [9l.Eo |, a uasidvnamical [ll| and partial dynamical symmetries [l3]) and 
their empirical evidence |l3l . [14L the study of finite-N effects 15|, [la, [l7|, [l8| and of the 



corresponding scaling behavior |l9l. l20l. l2ll. |22[|. The latter plays an important role in nuclei 



since the nucleus is composed of a finite number of particles, A^, and phase transitions, i.e. 
discontinuities in some quantities, are defined only in the limit A^ — t- oo. Also, the concept of 
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QPTs has been enlarged to include excited states quantum phase transitions (ESQPT) 
that is qualitative changes in the structure of a physical system as a function of excitation 
energy. Quantum phase transitions in the IBM are particularly interesting since this model 
has an algebraic structure, U(6), which is rather complex, giving rise within its parameter 
space to both first and second order transitions. On the other side, because of its algebraic 
structure, one can work out the "phases" in explicit analytic form, since these correspond 
to dynamic symmetries of the Hamiltonian. This situation is summarized in several review 



papers [2J, 125|, 126| , where a complete list of references is given. 

In this article we address the problem of how QPTs are affected by the presence of 
fermions in addition to bosons (QPTs in Bose- Fermi systems, odd-even nuclei). Nuclei 
offer a unique opportunity to study this problem since one has a model, the Interacting 
Boson- Fermion Model (IBFM) [i^l, in terms of correlated pairs with L = 0, 2 (s,d bosons) 
and unpaired particles with angular momentum j {j fermions), where this problem can 
be addressed in explicit form. (The problem is also of interest to the condensed matter 
physics community where one is interested in the motion of fermionic impurities in a bath 
of bosons [28] .) Studies of QPTs in odd-even nuclei were implicitly initiated years ago by 



Scholten [29|. Several explicit studies have recently been made by Alonso et al. |30l. l3ll. l32 



and by Boyiikata et al. [33[, who also have suggested a simple form of the IBFM Hamiltonian 
particularly well suited to study QPTs in odd-even nuclei because of its supersymmetric 
properties. Here we further expand on these studies and present the general theory of a 
single fermion with angular momentum j interacting with a system of s, d bosons (pairs 
with angular momentum L = 0,2) and study QPTs in this system both classically and 
quantum mechanically. We consider specifically the case of a particle with j = 11/2. After 
introducing the model Hamiltonian in Section |2l we present in Section [3] a classical analysis 
with novel results on single particle levels in a deformed field with both /3 and 7 deformation 
of interest not only to the algebraic description but also to its geometric counterpart. In 
Section m we present a quantal analysis of the same problem, introduce correlation diagrams 
for Bose-Fermi systems and study the classical-quantal correspondence. 

QPTs in odd-even nuclei would be a mere academic (and very complex) exercise were not 
for the fact that there is experimental evidence for the occurrence of these phase transitions 
in the odd-proton nuclei eiPm, gsEu and 65Tb. In the last section (Section E]), we discuss 
this evidence and compare the data with realistic IBFM calculations. The evidence is 
particularly clear in the unique parity negative parity levels originating from the /iii/2 level 
{i = 5,j = 11/2) and we therefore study these states. The large value of j has also the 
advantage that it provides a way to unravel the many complications of QPTs in Bose-Fermi 
systems which are not present for small values of j, j = 1/2 being trivial and j = 3/2 being 
analytically solvable. 



2. Model Hamiltonian 

We consider here the model Hamiltonian of a system of bosons with angular momen- 
tum L = 0,2 {s,d bosons) coupled to a fermion with angular momentum, j, as exemplified 
in the IBFM l27|. 
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Here the subscripts B, F and BF refer to the boson, fermion, and boson- fermion parts of 
the full Hamiltonian and 
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In these formulas, dots ■ denote scalar products, symbols x denote tensor products and 
: denotes normal ordering. Also, s'^,djj^ {s,d^) (/i = 0,±1,±2) denote creation (anni- 
hilation) operators for s,d bosons and aj_m(aj,m) {fn = ±1, if, •••5 ±j) creation (anni- 
hilation) operators for fermions with angular momentum j. The adjoint operators are 
d =(—Yd d- =l—V~"^n- 

The quantum phase transitions of the boson part of the Hamiltonian ([2]) are very well 
known |3J|. There are three "phases", characterized by their symmetry, U(5), SU(3) and 
SO (6), and two control parameters, ^ and x- We consider here the range of parameters, 
0<^<1, — ^<X<0. As^ changes from to 1 the system undergoes a quantum phase 

transition. No phase transition occurs when x changes from to — The phase transition 
as a function of ^ is first order from U(5) to SU(3) (spherical to axially deformed phases) and 
second order from U(5) to SO (6) (spherical to 7- unstable deformed phases). The situation 
is summarized in Fig. [1], shown here for sake of later discussion. 

In this paper we are interested in what happens when we vary the control parameters in 
the Bose- Fermi coupling, A, F, A, of Eq. ([3]). This problem is similar to (but more complex 
than) that of a bosonic system in an external field already discussed years ago by Landau (35| 
where one is interested in what happens as one varies the strength of the field, and to that 
of a fermion in a bath of harmonic oscillator bosons recently discussed in condensed matter 



physics [28 



S0(6) 




U(5) SU(3) 

Figure 1: Phase diagram of a system of s, d bosons. The three phases are denoted by their symmetry, U(5), 
SU(3) and SO (6). A hne of first order transitions ends in a point of second order transition in-between the 
U(5) and S0(6) hmits. The spinodal, critical and antispinodal points are denoted by respectively. 



3. Classical analysis 

Since the system we are considering here is that of an ensemble of bosons coupled to a 
single fermion, it is convenient to analyze the situation in terms of the motion of the single 
particle in the external field generated by the bosons This is done by introducing a 
(number projected) boson condensate and evaluating the expectation value of the Hamilto- 



nian H in the condensate [36|, |37| . This produces the fermion single particle Hamiltonian 
in the boson field, H. Diagonalization of H gives the single particle energies. The total 
energy is the sum of the single particle energy and of the boson energy. Minimization of 
this total energy with respect to the boson classical variables gives the equilibrium values 
(classical order parameters). The behavior of the order parameter (s) as a function of the 
control parameter(s) determines the nature and order of the phase transition. 

3.1. Expectation value of Hb in the boson condensate 

We begin by considering the boson condensate with good particle number N (3. 38, 39 



\N;/3,j) = -7=_mP,l)]"'\0) (5) 
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The expectation value of the boson Hamiltonian Hb of Eq. (|2]) in this condensate is given 
by (3 
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5.^. Expectation value of Hp and Vbf in the boson condensate 

By integrating out the boson degrees of freedom, i.e. by taking the expectation value of 
Hp and Vbp in the boson condensate, one obtains the fermion Hamiltonian 
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The matrix gmi,m2 is a real, symmetric matrix, with explicit form given by 

9m,,m2{N; /3, 7) = <X(iV; /3, 7) + gZi^^i^'^ 7) + ^7™^^,(iV; /3, 7) 

and 
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(The values for mi < ma can be obtained by noting that the matrix is symmetric. The 
symbol (.. | ..) in Eqs. (1121) and (IT^ denotes Clebsch-Gordan coefficients and the curly 
bracket in (fH]) is a 6j-symbol.) Also, in Eqs.(ll)-(13) we have isolated the dependence on 
N and on the parameters A, F, A, from that on the intrinsic variables /3 and 7 by introducing 
the matrices g^^ON ^ -quad ^ ~exc _ rpj^^ analysis of this section becomes then independent of 
the values of the couplings parameters and of N. 

3. 3. Diagonalization of the g-matrix 

The basis for the diagonalization of g is the fermion single-particle basis 
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Because of the structure of g, which couples basis states differing by ±2 in m, the diagonal- 
ization of g splits into two (doubly degenerate) pieces with m = j,j — 2,j — 4, — (j — 1) 
and similarly with m — )■ — m. The dimension of the basis is thus J + |- (In the case we 
discuss here, j = 11/2, the matrices are 6x6. When j = 3/2 the matrices are 2x2 and the 
results of the diagonalization can be obtained in explicit analytic form. When j = 1/2, the 
matrix is 1 x 1 and only g^'^^ and the Xq term in g^^^ contribute trivially). 

The diagonahzation of the matrix ^mi,m2(/3,7) yields the single particle eigenvalues 



^i(/3,7;x) 



1,2, 



1 



(16) 
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and eigenfunctions 

m 

These are the single particle levels in the deformed (3 and 7 field generated by the bosons. 
They depend on the boson parameter x which distinguishes different types of boson fields 
(x = 0, S0(6) symmetry; x = SU(3) symmetry). 

The contribution to the classical energy coming from the Bose-Fermi interaction is the 
single particle energy 

e,(iV; /3, 7; x; ^, T, A) = iV [^ef '^^(Z?, 7; x) - Tef 7; X) + A^f'^l/?, 7; x)] • (18) 

This contribution is proportional to the number of bosons iV as it should since the Bose- 
Fermi interaction, Vbf-, is linear in the generators of U(6). 

3.3.1. Results for the single particle eigenvalues 

(a) The monopole term. This term is diagonal with eigenvalues which do not depend on 
the index i and on the parameter x 

^?'°''(/J.7;x) = (^) . (19) 

(b) The quadrupole term. The eigenvalues of the matrix must be obtained nu- 
merically. A program has been written in Mathematica to diagonalize the matrix. Using 
this program, we have first investigated the symmetries of the eigenvalues. Fig. [2] shows the 
eigenvalues ef'^^^{(3,'j; x) as a function of 7 in the domain < 7 < 27r for x = — ^ and 
/3 = a/2 (top part) and for x = 0, /3 = 1 (bottom part). (We have chosen these values of 
f3 because they are the equilibrium values, /^e.B, of the boson energy functional, Eb^P,!), 
Eq. ([H]), for SU(3), x = and S0(6), x = 0). The various states are labelled by the 
index i = 1 , . . . , 6 in order of increasing energy. This figure shows the remarkable result that 
the single particle energies generated by the quadrupole Bose-Fermi interaction are periodic 
in 7 with period ^ and are symmetric under refiections around |. This is in spite of the 
fact that g contains terms cos 7, sin 7, cos 27, sin 27. We can therefore restrict our study to 
< 7 < f . (The periodicity ^ can be proven by constructing the characteristic polynomial 
that diagonalizes the matrix and showing that it is a function of cos 87). From Fig. |2]one 
can see that the 7 dependence of the single particle levels is different for x = than for 
X = 0. 

In Figs. [3] and m we show the eigenvalues ef^"^^ {(3,'j; x) as a function of /3 (at 7 = 0°), 
in the interval —2 < f3 < 2 (top part), and as a function of 7 in the interval 0° < 7 < 60° at 
l3 = \/2 for X = and at /3 = 1 for x = (bottom part). The bottom part is identical 
to Fig. |2]but restricted to 0° < 7 < 60°. We see that when x = we have an additional 
symmetry 

^r^'^(/3,7;0) = -er^^(-/3,7;0). (20) 
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Figure 2: Dependence on 7 of the eigenvalues sf^^^ (/?,7;x) (* = l,---7 6) of the matrix for 
X = — /3 = \/2 (top part) and x = 0, /3 = 1 (bottom part), labelled in order of increasing energy. 



This symmetry is lost when x 7^ 0. By combining the symmetries in 7 with those in /3 we 
find that the single particle energy eigenvalues 7; x) ^-re invariant, when x = 0, 

under the transformation 

^QC/AD^_^^ 7 = 60°; 0) = 7 = 0°; 0) . (21) 

This symmetry is lost when x 7^ 0- order to recover it, we need to simultaneously change 
the boson field from x < (prolate) to x > (oblate) when we go from /3 > to /3 < 0. 
The single particle energies have then the property 

^^^(-/3, 7 = 60°; -x) = 7 = 0°; x) • (22) 

For the purpose of this article, where we are interested in the effect of a single fermion on a 
bosonic system with a definite value of x < (prolate boson field) it is sufficient to restrict 
all further analysis to the domain /3 > 0, 0° < 7 < 60°. 
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Figure 3: Eigenvalues ef (/5,7;x) foi' ^ particle with angular momentum j = 11/2 in a quadrupole /3,7 
field with X — (prolate field) as a function of /3 for 7 = 0° (top part) and as a function of 7 for /? — \f2 
(bottom part). The eigenvalues are labelled by the projection of the angular momentum on the 3 axis at 
7 = 0° and on the 2 axis at 7 = 60°. 



We also note that the matrix at 7 = 0° is diagonal. Its eigenvalues are given in 

terms of the projection of the angular momentum on the intrinsic axis 3. The eigenvalues 
at 7 = 60° are instead given in terms of the projection K2 on the intrinsic axis 2. There is 
a continuous change from K^^ to K2 with the correspondence shown in Figs. [3] and HI The 
single particle levels in a /3, 7 field were studied years ago by Meyer-ter-Vehn [40| within the 
framework of the collective model. The 7 dependence of our results for x = is identical to 
that of 40| . On the other side the /3 dependence differs for large /3 due to the denominators 
Yqig2 in our ^-matrix and to the fact that the variable (3 of the collective model, I^bmi is related 
to that of the Interacting Boson Model, Pibm, by a scale transformation, Pbm = c/3ibm with 
c ~ 0.2 in the mass A = 150 region [2i]. We also note that, if the variable 7 is frozen to 0°, 
one obtains the single particle levels in a deformed /3 field with axial symmetry. Those were 



studied years ago by Nilsson within the framework of the collective model (4l|. Our results 



9 



QUAD 





0.8 








r=o° 






Ji2___- 


0.4 






5/2 






7/2 


7/2 






^_J/2_ 


~9/2~~~ 






— — — ^ 










-2 


1 


2 



QUAD 

/■ 

0.8 r 




Figure 4: Eigenvalues iPillX) ^oi a particle with angular momentum j = 11/2 in a quadrupole /3,7 

field with x = (7-unstable field) as function of j3 for 7 = 0° (top part) and as a function of 7 for /3 = 1 
(bottom part). Labelling as in Fig. [31 



for X = ^re similar to those of the Nilsson model except for the denominators j^^- 

(c) The exchange term. The eigenvalues of the matrix gmfm2 obtained numer- 
ically. We first investigate the symmetries of the eigenvalues. Fig. |5] shows the eigenvalues 
ef'^'"(/3, 7; x) as a function of 7 in the domain < 7 < 27r for % = — = -\/2 and 
X = 0, /3 = 1. This figure shows that the contributions to the single particle energies coming 
from the exchange interaction are periodic in 7 with period |, that is half of the period 
of the quadrupole interaction, and are symmetric under reflection around |, i.e. they are 
a function of cos 67. Furthermore, the exchange term is independent of x and is an even 
function of (3. 

(d) Combination of quadrupole and exchange. We consider here the eigenvalues, 
^QC/AD+Exc^^^ ^. ^£ matrix 



C^r^''^(A7;x) = -C;:.?(/3,7;x) +0;) ~9^f,ZMr,x) (23) 



10 



EXC 
/■ 




-0.4- 

60° 120° 180° 240° 300° 360°'/ 

Figure 5: Dependence on 7 of the eigenvalues ef-^'~" {(3,'-f;x) (* = l7---i6) of the matrix g^^^^^ for x — 
— ^,13 = \/2 (top part) and x = Oi/^ = 1 (bottom part), labehed in order of increasing energy. 



and, for sake of display, we take (^) = 3. The eigenvalues of this matrix are shown in Figs. |6] 
and [71 We see from these figures that the ordering of the single-particle eigenvalues at a 
given f3, 7 when both quadrupole and exchange interactions are present is rather complex. 
Our results are similar to those of the Nilsson model plus BCS 27|. The similarity is no 



longer so obvious as in the case in which the coupling is purely quadrupole, but it can be 
seen by solving the BCS equations in the deformed /3, 7 field as done by Meyer-ter-Vehn 4G . 



The formalism is flexible enough that one can also investigate particle- hole conjuga- 
tion. Within the framework of IBFM, particle-hole conjugation is the transformation F — )■ 
— F, A — )• A (see the following Eq. fl38p ). When A = 0, the effect of this transformation is to 
reverse the order of the single particle levels in the previous Figs. |3] and HJ as one can simply 
see from the form of g^^^^[[3,'~^; x) iii Eq. f lT2|) . When A 7^ the transformation leads to 
a different ordering of the single particle levels. We do not dwell on this point further, but 
note instead that a special circumstance occurs when the variable 7 is frozen to 0° (prolate 
axial symmetry). As remarked above, the quadrupole Bose- Fermi interaction is diagonal at 
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Figure 6: Eigenvalues ef (/?i7;x) for ^ particle with j = 11/2 in an external /3,7;x — ~ 2 ^^^'^ 

with both quadrupole and exchange interaction with strengths A/F = 3. Labelling as in Fig. [31 Figure 
constructed by using F = -0.0958, A = -0.287. 



7 = 0°. It turns out that also the exchange interaction is diagonal. The eigenvalues of the 



5f-matrix for the combined quadrupole plus exchange interaction are given by |36 





where 



-NK { [ ) (2j + l)P^[^K^ - j{j + 1))]^ 1^ (24) 



P, = [{2j - l)j(2j + l)(j + l)(2j + 3)]-^/' . (25) 



The eigenvalues can be labelled by the projection of the angular momentum on the intrinsic 
axis 3, K3 = K = ^, |, j and they are doubly degenerate. 

The eigenvalues of the g-mathx for the combined quadrupole plus exchange interaction 
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Figure 7: Same as Fig. [S] but with x = 0. 
are also diagonal at 7 = 60° (oblate axial symmetry), and are given by 

-NA I (2j + l)Pn^Kl-j{j + l))]^! , (26) 

where K2 is the projection of the angular momentum on the intrinsic axis 2. They satisfy 
the relation 

ookM = ^K,^K,{-^) . (27) 

3.4- Equilibrium values: classical order parameters 

Having determined the single particle energies, e,, Eq. f|T8|) . we can construct the total 
energy surface 

Ei{N; 13, 7; e, x; A, r. A) = EsiN; /5, 7; x) + e, + e,(iV; /5, 7; x; A, T, A) , (28) 
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Figure 8: Equilibrium values, j3e,B, as a function of the control parameter, of the bosonic part of the energy 
functional, Eb, for the transition U(5)-SU(3), X = left panel, and U(5)-SO(6), x = 0, right panel. 

The spinodal, critical and antispinodal points are denoted by C'^, ^**, respectively. Figure constructed for 
iV = 10 bosons. 



where Eb, Eq. ([7]), and hence Ei depend also on the scale parameter Sq. In the traditional 
approach to odd-even nuclei, the equilibrium values /3e,7e are obtained by minimizing the 
boson energy Eb (or core energy in the collective model). This provides a static deformation 
and the single particle energies are evaluated in this static deformation. This is a good 
approximation when the deformation is large. In the transitional region, this is no longer 
the case. A better approximation is to minimize the total energy surface including the 
contribution of the fermions (the odd-particle). Although the fermion contribution is of 
order it will be shown that it has a dramatic effect on the phase transition, modifying 
the location of the critical point. 

The equilibrium values /3e and 7e (the classical order parameters) for the combined system 
are obtained by minimizing Ei with respect to /3 and 7, i.e. imposing the conditions 

dEi dEi , , 

Minimization of Ei is, in general, no simple matter. The boson part, Eb-, has a minimum 
at (3e,B{0 and at 7e,B = 0°, for any value of C, and x < 0, except for x = where there is 
no minimum in 7 (7- unstable situation). The boson equilibrium values Pe,B{0 have been 
calculated by several authors [l6j and are shown in Fig. [8] as a function of the control 
parameter ^ for two values of x = and x = 0, for later comparison with those in 

odd-even nuclei (Bose- Fermi system). 

The single particle eigenvalues have a rather complex structure (Figs. |3]and|l]), espe- 
cially when there is an exchange interaction (Figs. Eland [7]). When doing the minimization, 
an important question is on what space to minimize. In view of the periodicity in 7, it is 
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sufficient to minimize Ei in the sextant 0° < 7 < 60°. Also, since the purpose of this article 
is to understand the extent to which the presence of an odd fermion modifies the phase 
transition, we seek solutions in the neighbourhood of the boson equilibrium values and thus 
we do a constrained minimization with /3 > 0, in practice < /3 < 2. 

To fully investigate QPTs in IBFM, one should do a study as a function of all con- 
trol parameters ^, x (bosons) and y4,r,A (Bose- Fermi). In this paper we study two phase 
transitions: (i) U(5) to SU(3), spherical to axially deformed, X = ^"^^ ip) U(5) to 
SO (6), spherical to 7- unstable, x = 0, both induced by a change in < ,^ < 1. These 
phase transitions have been investigated extensively in the purely bosonic case and are of 
practical importance. We also set A = 0, since the monopole term only renormalizes the hd 
term in the boson Hamiltonian, Hb-, and its effects can be easily seen. We are left with two 
additional control parameters, F and A. In order to isolate the dependence on the fermion 
angular momentum j, we write 



r = Fog 



Qn = {j II y^^^ II j) = V v/2j + 1 (j, 1/2; 2, | j, 1/2) (30) 

where the double-bar in Eq. (1301) denotes a reduced matrix element. To further restrict the 
study, we take Fq proportional to ^ 

To = . (31) 

The entire study presented here is thus in terms of a single control parameter ^, plus the 
value of A which, for convenience, can also be written as 

A = Aeo (32) 

making all terms in H then proportional to the scale factor Eq. The parameterization 



of Eq. ([21]), for the quadrupole term, was suggested by Alonso et al. [30[ and it has the 
advantage that the combined Hamiltonian can be rewritten, apart from an overall constant, 
as 



H = 5o|(i-0(^rf + %)-^(0^ + ^F)-(g^ + gF^ 

+A v/2jn : [{S X a,)(^) x {d x aj)(^)] ^ : | , (33) 



where 



Qf = Qjj Qj , (34) 



and Qj = {Oj X dj)^'^\ If one multiplies Qjj by —y/n, the corresponding Hamiltonian for 
j = 3/2 and x = A = 0, has a Bose- Fermi symmetry SO (6) SU(4) D Spin(6) [i^], i.e., the 
Hamiltonian is invariant under interchange of the boson operators ri^, with the fermion 
operators fij, qj. 
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Figure 9: Equilibrium values, /3e,i (top part) and ■je.i (bottom part) as a function of the control parameter 
^ in the U(5)-SU(3) transition (x — —^), A = and TV — 10. States are labelled by the index i — 1, 6. 



3.4.1. Results for A = 

The equilibrium values for the U(5)-SU(3) transition are shown in Fig. inland, for the 
U(5)-SO(6) transition, in Fig. (TDl From these figures one can see the effects of the odd 
particle on the phase transition. In the U(5)-SU(3) case, the phase transition is washed out 
for states 1,2,3 and enhanced for states 4,5,6. The critical point is approximately at the same 
location as for the purely bosonic case (^^ ~0.51) for state 4 (^| ~0.50), but it is moved 
to larger values for state 5 ~0.53) and state 6 (^g ~0.58). The values of ■je below and 
around the critical point are no longer zero. After the critical point, all states become those 
of a single particle in an axially deformed field (7e = 0°) with equal deformation I3e,i = (3e,B- 
States in this region can be labelled by the projection K = 1/2, 3/2, 5/2, 7/2, 9/2, 11/2 of 
the angular momentum on the intrinsic 3 axis, corresponding to states i = 1, 6. A different 
situation occurs for the U(5)-SO(6) transition. Although the equilibrium values, /3e,j, follow 
the same behavior as for the U(5)-SU(3) transition, the equilibrium values 7e,j do not. 
States 1 and 6 are oblate {■je = 60°), states 2,3,4 are prolate (7e = 0°) and state 5 is triaxial 
(7e = 30°). For this phase transition, the equilibrium value of 7 is dictated by that of the 
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Figure 10: Equilibrium values, Pe,i (top part) and ^e.i (bottom part) as a function of the control parameter 
^ in the U(5)-SO(6) transition (x'= 0), A = and N = 10. Labelling as in Fig. El 

single particle states, Cj, since the boson part is 7 independent. The effect of the odd-particle 
is more dramatic here than in the case of the U(5)-SU(3) transition. 

3.4.2. Results for 

The equilibrium values when A = —0.287 for the U(5)-SU(3) transition are shown in 
Fig. [TT]and for the U(5)-SO(6) transition in Fig. [121 We see that the modifications induced 
by the presence of a fermion to the phase transition are more dramatic when the exchange 
interaction is added, especially in the location of the critical point The state 6 moves 
considerably to the right of In summary, QPTs in the presence of an odd fermion have 
a different behavior as a function of the control parameter ^ than in the purely bosonic 
system. The odd fermion acts as a catalyst for some states and as a retardative for others. 

3.5. Single particle energies as a function of the control parameter ^ 

Once the values of /3e,j and , 7e,i have been obtained, one can evaluate the single particle 
energies ei{N; /3e,i, 7e,i) as a function of ^. Throughout this section we set the scale parameter 
£0 = 1 and = 10. 
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Figure 11: Same as Fig.[ibut with A = -0.287. 



3.5.1. Results for A = 

Results for X = X = ^ind A = are shown in Fig. [131 The discontinuities reflect 

the phase transition. 

3.5.2. Results /or A ^ 

Results for X = — X = and A = —0.287 are shown in Fig. [TH The single-particle 
energies are more spread in this case than when A = and have a different ordering. 

3. 6. Total energies as a function of ^ 

It is of interest to show also the results for the total energies Ei{N; I3e,i, 7e,j) as a function 
of ^. (The total energies depend also on x and A). Throughout this section we set the scale 
parameter Eq = 1 and = 10. 

3.6.1. Results for A = 

The effects of the odd particle in the total ground state energies are small, of order 1/A^, 
as shown in Fig. [151 
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Figure 12: Same as Fig.[Tn]but with A = -0.287. 



3.6.2. Results for A ^ 

The effects of the odd particle are also small here, as shown in Fig [161 for ^ = —0.287, 
although larger than in the case of A = since the reordering of the single particle states 
produces a lowering of some of them. 



4. Quanta! analysis 

The crucial aspect in the study of quantum phase transitions is obviously the quantal 
analysis. This analysis must match the classical analysis of the previous section, valid in 
the limit N ^ oo. The quantal analysis for IBFM is done by diagonalizing numerically the 
Hamiltonian 

H = Hb + Hf + Vbf (35) 

for finite N using the program ODDA ^]. We take in this study N = 10. Although this 
is a relatively small value, previous studies in the purely bosonic case have shown that the 



salient features of QPT are already apparent [16] at = 10. 
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Figure 13: Minimum single particle energies, {ei)^^^ = —NT {si)^^^ as a function of the control parameter, 
^, in the U(5)-SU(3) transition (x = (top part) and in the U(5)-SO(6) transition (x = 0) (bottom 

part), both with A = 0. 

The program ODDA uses a semi-microscopic version of IBFM in which (for a single 
fermion with angular momentum j) 



and 



Vbf 



^3 nj 



■irMON I -i^QUAD I -i^EXC 

Vbf "t ^bf "t Vbf 



(36) 



V, 



MON 
BF 

QUAD 
BF 



As fid hj 



r, (m] - v]) Qj, ■ q, + q, ■ Q> 



(37) 



with Qjj defined in Eq. fISU]) . In this semi-microscopic version, the Bose- Fermi interaction is 
given in terms of the BCS occupation probabilities, Uj and vj, with u'j + Vj = 1. Comparing 
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Figure 14: Same as Fig-Obut with A = -0.287. 

with the model Hamiltonian ([3]) of Section [21 we have the relationships 

A = A, 

r = r, 2 («| - Q,,- 

A = -k,%^u]v]Q%l(23\V). (38) 

J^.l. The transition from spherical to axially deformed (U(5) to SU(3)) 

We study this transition by using a slightly modified form of the boson Hamiltonian 



H 



U{5)-SU{3) 



^0 



il-OCiiU{5))-^^C,{SU{3)) 



(39) 



with X = III (139]) a term L ■ L has been added to the Hamiltonian of Eq. ([2]), where 

L is the boson angular momentum operator. The combination 2[Q ■ Q + ■ L] is then 
the quadratic Casimir operator of SU(3). Also we denote by Ci{g) and C*2(fi') Casimir 
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Figure 15: Minimum total energies, Ei{N; /?, 7; x) as functions of the control parameter ^, in the U(5)-SU(3) 
transition (x = (top part) and in the U(5)-SO(6) transition (x = 0) (bottom part), both with A = 0. 

operators of the first and the second order of the algebra g. The reason for using this boson 
Hamiltonian is that we want to isolate the intrinsic from the rotational part of the spectrum, 
in order to make the comparison between classical and quantal analysis straightforward. In 
the Bose-Fermi interaction, we set Ag = for the same reason as in the previous section. 

4-1.1. Correlation diagrams 

We first study the nature of the spectra in the presence of only a quadrupole interaction, 
V^p^^ , and set, 

r. = ^eo. (40) 

In the parameterization (l37|) we need to specify the values of the occupation probabilities, 
u|. The correlation diagrams for = 0,f| = 1 (particle-hke spectra) and m| = = 
(hole-like spectra) are shown in Fig. [T71 Correlation diagrams describe how the energy levels 
evolve from one phase to the other as a function of the control parameter ^. The phases 
are defined here by their symmetries, U(5) spherical, SU(3) axially deformed, and S0(6) 
7-unstable. Fig. [T7] shows how rotational bands emerge from the spherical basis. It also 
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Figure 16: Same as in Fig.HUbut for A = -0.287. 



displays in a clear fashion particle-hole conjugation, that is the transformation [u < — > v). 

With the values of the BCS coefficients used in Fig.[T71 the exchange interaction vanishes. 
To study the effect of the exchange interaction, we construct the correlation diagram as a 
function of f|. Fig. [13 This is the same as varying A, since m| = 1 — f | and thus 

A = -A, 8v^(l - v])v]Q'^^/i2j + 1) . (41) 

As Vj changes from to 1, we span the entire set of spectra from hole- like to particle- like. 
Particularly interesting in this figure are the spectra at half-filling, Vj = 0.5. At this point, 
the quadrupole interaction vanishes, and the spectrum is given by the exchange interac- 
tion. It should be noted, however, that the diagram shown in Fig. [18] is an unconventional 
correlation diagram, since by varying Vj we vary both the exchange and the quadrupole 
interaction, Eq. (137|) . 

4.1.2. Classical- quantal correspondence 

In addition to providing a description of the evolution of the quantal levels from one 
symmetry (phase) to another, correlation diagrams are also useful to test the classical- 
quantal correspondence. In Section [3l the single-particle energies for a fermion with angular 
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Figure 17: Correlation diagram for a j ~ 11/2 particle coupled to a system of {s,d) bosons undergoing a 
U(5)-SU(3) phase transition. Top part v'^ — 1, bottom part = 0. The interaction is purely quadrupole. 

momentum j in the presence of the boson condensate of Eq. ([6]) were calculated classically. 
For 7 = 0° they were given by Eqs. fl24p and fl25p . as a function of f3 and for any x- When 
the bosons have SU(3) symmetry, X = iii the limit — ?■ oo, the equilibrium 

deformation is /3e = \/2. Inserting these values in Eq. fl21]) we obtain the classical result 

-NA^ (2j + 1) - jU + (42) 

An analytic quantal calculation of the single-particle energies for bosons with SU(3) symme- 
try in the limit of large N can also be done. The bosonic ground state wave function, with 
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Figure 18: Correlation diagram for a. j — 11/2 particle coupled to a system of (s, d) bosons with quadrupole, 
Vgp^^,X — and exchange, V^p'" , interaction parametrized as in Eq. (I57|) . with As/F^ = 3, as a 

function of the filling probability = v'j of the single particle state. States are labelled by the total angular 
momentum J. 

good SU(3) quantum numbers, can be written as | [A^] , {2N, 0) ,Kc = 0, L, M) [7]. States for 
a fermion with angular momentum j coupled to the bosonic ground state can be written 
as I [A^] , (2A^, 0), = 0; j, Kj; K = Kj, J, M). A long but straightforward calculation of the 
matrix elements of the interaction Vbf = V^p^^ + V^p'^, with V^p^^ and V^p'^ as in 
Eq. dSZ]), gives H 



hm m , (2Ar, 0) , 0; j, K,- K,J,M\ Vbf \ [N] , {2N, 0) , 0; j, Kf, K, J, M) = Xf^'^ . (43) 



N- 



Here the parameters T and A in the expression for Xf^'\ Eq. m, 

are related to the 

parameters of Vbf, Eq- (|37|) . by means of Eqs. (138|) and (HOl) . Specifically, T = 2Ts{l — 
2Vj)Qjj, Tg = ^ Eq and A is given in Eq. (HT]) . It should also be noted that in the limit 
N ^ oo, the expectation value of Vbf is independent of the the total angular momentum J. 
Both the explicit formula Eq. fj43|) and the independence on J can be checked numerically. 



In Fig.[T7]top, obtained with T = — ^ (2£:o) Qjj and A = 0, one can see that at ^ = 1, SU(3) 
symmetry, all rotational levels built on the intrinsic state K, converge to a single value, as 
seen on the right-hand side of the figure, and that the energies of the intrinsic states can be 
described to a good approximation by Eq. f H2]) . Since the figure represents the results for 
A^ = 10, one can see that small deviations occur, especially for the band with K = 1/2, but 
these deviations are relatively small and the large A^ limit appears to be already reached. 
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The same conclusion applies to Fig. [T7]bottom obtained with T = +^ (S^o) Qjj, and A = 0. 
Flipping the sign of F (particle-hole conjugation) reverses the ordering of the K states but 
the energies are still given by Eq. ( H2|) . (The bottom part of Fig. [T7]is also consistent with 
the top part of Fig. [3l once the relationship F = 2TsQjj is taken into account). 

The correspondence between classical and quantal calculation is exact at ^ = 1, SU(3) 
symmetry, and — ?■ oo. As one moves away from ^ = 1, or for finite A^, it becomes 
approximate. Nonetheless, rotational bands can still be identified for C,c ^ ^ ^ 1- Below the 
critical value, the characterization of states by a A' quantum number is no longer possible. 
For this region, < ^ < a straightforward comparison between classical and quantal 
calculation is no longer possible and one must resort to a numerical calculation both classical. 
Section [31 and quantum-mechanical, Section HI 



4.1.3. Ground state energy 

The ground state energy is a key indicator of phase transitions. The ground state energy, 
£"0, its first, and second, ^p^, derivatives with respect to ^ are shown in Fig. [19] (particle- 
like, f| = 1). (The energy Eq in this figure is in units of the scale factor taken to be 
£0 = 1 and A^ = 10). For hole-like spectra, Vj = 0, we have a similar behavior, as well 
as in the case in which A 7^ 0. The quantal result for Eq should be compared with the 
classical result, {Eo)^iri, shown in the top part of Fig. [151 where (-Eo)min is the energy of 
the lowest state at the minimum. Within 1/A^ corrections, the two results agree (classical- 



quantal correspondence). Fig. [THjis similar to Fig. (7.3) of j2l|] for the purely bosonic case. 
Although the value of A^ used here is too small to distinguish between first and second order 
transition, with discontinuities in and respectively, nonetheless precursors of the 
QPT are clearly seen at ^ = — 0.5. 

4.1.4. Quantal order parameters 

The equilibrium deformations Pe,i,1e,i are the classical order parameters. As quantal 
order parameters we consider here only the expectation value of in the states i = 1, 6, 
(first quantal order parameter), z/}^"* 

(1) ^ {^P^\na\i^^) , . 

' N ' ^ ' 

This is shown in Fig. [20] top part. Since the order parameters z/|^'' are related to the square 
of the classical order parameters /3e,i, this figure is related to Fig. [9] to which it corresponds 

(1) di/^^'^ 

in the limit A^ — t- 00. The derivative of in the ground state, is also shown in 
Fig. [20] bottom part. (This quantity diverges when A^ — t- 00). From this figure one sees 
clearly that the transition is made sharper by the presence of the fermion for some states, 
11/2, 9/2, 7/2, while is made smoother for others, 5/2, 3/2, 1/2, a result already seen in the 
classical analysis. 
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Figure 19: The ground state energy, Eq, top, its first derivative, center, and its second derivative, 
bottom, for a j = 11/2 particle coupled to a system of {s,d) bosons undergoing a U(5)-SU(3) transition. 



4-2. The transition from spherical to 'y— unstable (U(5)-SO(6)) 

We study this transition by using the standard form of the transitional boson Hamiltonian 

(45) 

with X = 0. We set = as in the previous section and Vg = -^^o as in Eq. fHO|) . 

4-2.1. Correlation diagram 

The correlation diagram for this phase transition is shown in Fig. In this case, there 
is no clear separation between intrinsic and rotational motion, since Q^^^ ■ Q^^^ is not a 
Casimir operator of SO (6), and thus the states at ^ = 1 do not coalesce into a single point 
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Figure 20: The quantal order parameters, , Eq. (|44p . as a function of the control parameter, ^ (top 
part). The expectation values of the number of d bosons, {fid), for the lowest states of a given total angular 
momentum, J — 11/2,..., 1/2 for a j = 11/2 particle coupled to a system of {s,d) bosons undergoing a 

g^(l) 

U(5)-SU(3) transition. Particle-like spectra, v"^ — I, and no exchange interaction. The derivative as a 
function of ^ (bottom part). 

as in Fig. [T71 (We have also done a study in which C2{SO{6)) is used instead of Q^"" ■ Q^^^- 
These two operators are related by (72(5*0(6)) = Q^^^ ■ Q^^^ + C*2(>S'0(5)). However, using 
the Casimir operator, makes the figure at the end point ^ = 1 very crowded, since all states 
arising from a given 5*0(6) representation collapse to zero energy, and, for this reason, we 
prefer to plot in Fig. [5T] the results with the Hamiltonian fHSl) . 

The effect of the exchange interaction is studied by constructing the correlation diagram 
as a function of f|. Fig. [221 A remarkable property of this diagram is the symmetry under 
particle-hole conjugation. The diagram is symmetric around half-filling. This is unlike the 
case of the U(5)-SU(3) transition. 
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Figure 21: Correlation diagram for & j — 11/2 particle coupled to a system of (s, d) bosons undergoing an 
U(5)-SO(6) transition. The interaction is purely quadrupole. There is no difference in this case between 
particle- like (w^ — 1) and hole-like (w^ = 0) spectra. 



4-2.2. Ground state energy 

The ground state energy, Eq, its first and second derivative are shown in Fig. [231 U(5)- 
S0(6) transition. This figure is similar to Fig. Uni U(5)-SU(3) transition, and as mentioned 
in Section I4.1.3[ it is not possible to distinguish whether the transition is first or second 
order. In order to do so, one must go to much larger values of N, as done in the purely 
bosonic case ^]. However, the fact that the quantity is smoother in Fig. [23] than in 
Fig. [m supports the conclusion that the ground state transition is first order for U(5)-SU(3) 
and second order for U(5)-SO(6), since the second derivative of Eq diverges (or not) for first 
(or second) order transitions in the limit N ^ oo. 

5. Experimental evidence 

Even-even nuclei in the mass region A~150 are known to experience a first order quantum 
phase transition (U(5)-SU(3)) at neutron number 90. The evidence for this statement comes 
from the analysis of (i) the excitation spectrum which displays a gap, A = -E'(O^) — £'(0^); 
(ii) the two neutron separation energies, S2n{N) = —[Eq{N + 1) — Eq{N)], proportional to 
the derivative of the ground state energy, Eq, with respect to the control parameter, ^, i.e. 
and (iii) the B(E2;0 — i- 2) values, proportional to the square of the order parameter 

{fid)^ [13, El, [13] • In odd-even nuclei, the first two quantities can be easily measured, while 
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Figure 22: Same as Fig. [18] but with % = 0. 



the last quantity, B{E2; Jg g — )■ J') is more difficult to measure due to the fragmentation of 
the B{E2) strength from the ground state to several states J'. In this paper, we therefore 
analyze the ffist two quantities and show that both display the features expected for a QPT, 
as described in Section |3] and Section HI 



5.1. Excitation spectrum 

Odd-proton nuclei in the region of A~150 offer an unique opportunity to study QPT 
in Bose-Fermi systems, because of the occurrence of the unique parity state /iii/2 near the 
Fermi surface. We therefore analyze the negative parity states of the odd-proton nuclei, 
6iPm86_92, 63EU86-92, and 65Tb86_92- A realistic calculation for these nuclei can be done 
within the framework of IBFM, with an Hamiltonian H = Hb + Hp + Vbf- The IBM 
Hamiltonian Hb is written as 



Hb = ehd+ J] cl^V2L + 1 {d^ x S)^^^ x (d x (tj 



L=0,2,4 



[d^ X rft)(2) X (dx sV'' + [s^ X d^y X (dxd 



(2) 



(2) 



(0) 



1(0) 



(rft X dt)W X (s X sY'^ + (st X xidxd 



A (0)' 



(0) 



(46) 
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Figure 23: Same as Fig. [12] but for the U(5)-SO(6) transition. 



In order to perform a calculation of the odd proton nuclei indicated above, we need the 
Hamiltonian parameters of the even-even nuclei 6oNd86-92, 62Sm86_92, 64Gd86-92- We take 
these parameters from previous studies. In these studies, first a calculation in the Proton- 
Neutron Interacting Boson Model (IBM-2) is done with parameters given in Table [1] The 
input parameters in IBM, e, cl(L = 0, 2, 4), t>2, fo, Eq. fH6|) are then calculated from these 



by a projection technique |43 



In the case of the unique parity configuration /iii/2 there is only one single particle energy, 
Ej, which we take as Ej = 0. The Bose-Fermi interaction ( 137|) is specified by the parameters 
As, Ts and and f|, three of which are independent. Using the semi-microscopic theory, 
one can extract the occupation probabilities The parameters used in the present analysis 
are given in Table El For gsEu and QiPm isotopes, they are the same parameters used by 
Scholten and Blasi |29| and by Scholten and Ozzello |4J|, respectively, in previous studies of 
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Table 1: Parameters of IBM-2 used in the calculation of the even-even nuclei i^*^-i^°Nd, "8-1523^^^^ 
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1.05 


1.97 


0.38 


155Tb 


-0.1 


1.26 


2.79 


0.40 


157Tb 


-0.1 


1.41 


2.87 


0.41 



Table 2: Strengths of the Bose-Fermi couplings in the odd-even nuclei ^"^"^-^^^Vm, ^^^^-i^seu, ^^i-^^'^Tb. 
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86 88 90 92 86 88 90 



Neutron number 



Figure 24: Comparison between calculated and experimental spectra of negative parity states in giPni. The 
lowest 11/2^ state is taken as zero of the energy. The parameters of the calculation are given in Tables [T] 
and [21 In the experimental spectra, taken from 45], uncertain assignments are indicated by open symbols. 



these isotopes. For the 65 Tb isotopes we scale the latter parameters. Our calculated spectra 
are shown in Figs. [2H [25l [26l where they are compared with the available experimental 



data (45| . One can see very clearly the phase transition occurring between neutron numbers 



88 and 90 both in the theoretical (left) and in the experimental (right) spectra. 

5.2. Two-neutron separation energies 

A complete analysis of two-neutron separation energies requires a calculation of both 
positive and negative parity states in QiPm, es Eu and esTb. For the Eu isotopes this 



calculation was done by Scholten [29| . The experimental two-neutron separation energies in 
6iPm, esEu and 65Tb |46| are shown in Fig. [271 

One can see clearly the occurrence of discontinuities in the behavior of the two-neutron 
separation energies, an indication of a QPT. In order to emphasize these discontinuities, we 
note that the two neutron separation energies are given by a smooth contribution linear in 
N plus the contribution of the deformation 0, p.74] 

S2n = -A2n-B2nN + St^ . (47) 

In Fig. [28] we show the deformation contribution, extracted from the data by subtracting the 
linear dependence, with = —15.185, —16.37, —17.672 MeV for Pm, Eu, Tb, respectively, 
and i?2n = 0.670 MeV. The deformation contribution can be easily calculated using IBFM. 
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84 86 88 90 92 84 86 88 90 



Neutron number 

Figure 25: Same as Fig. [24]but for gaEu. 




86 88 90 92 86 88 90 



Neutron number 

Figure 26: Same as Fig. [^Hlbut for 65Tb. 
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Figure 27: Experimental two-neutron separation energies, S{2n), for the even-neutron isotopes of eiPm, 
63EU and 65Tb, taken from 46 1. 



Fig. I29] shows this contribution for the hii/2 level, for QiPm, esEu and 65Tb. Figs. 1281 and 1291 
cannot be, in principle, directly compared since the experimental data are the separation 
energies of the actual ground state, while Fig. 129) shows the separation energies of the /iii/2 
level. However, in practice, differences between the two are of order and thus, to a good 
approximation. Fig. [2H] can be compared to Fig. [221 In particular, the onset of deformation 
is clearly seen both in Fig. and Fig. 



6. Conclusions 

In this paper, we have analyzed QPTs in Bose-Fermi systems, specifically the effect of 
one single-fermion with angular momentum j immersed in a bath of bosons with angular 
momentum L = 0,2 {s,d bosons). 

By doing a classical analysis (Section |3]) we have studied the single particle motion in a 
field with /3,7 deformation and determined the classical order parameters /3e,i,7e,i- We have 
shown that while the presence of the odd-fermion does not influence much the motion in 
the strongly deformed regions, the corrections being of order it does influence greatly 
the location of the critical point and the entire nature of the phase transition, washing out 
the transition for some states and enhancing it for others. This is a novel result which may 
have applications to other fields of physics outside of nuclear physics. 

By doing a quantal analysis (Section H]) we have investigated the change in level structure 
induced by the phase transition (correlation diagrams). Here we have obtained the result 
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Figure 28: The contribution of deformation, S{2n)def to the two-neutron separation energies, Eq. (|Tf)) . in 
6iPm, esEu and 65Tb. Experimental data taken from (46| . 
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Figure 29: The contribution of deformation, S{2n)def Eq. (j47p . calculated in IBFM with the parameters 
given in the text. 
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that while the presence of the odd-nucleon does not affect much the ground-state energy 
correction) it does affect greatly the level structure. This level structure is rather complex, 
especially in the neighborhood of the critical point. Wc have also used the correlation 
diagrams to investigate the classical-quantal correspondence and shown that when the bosons 
have SU(3) symmetry the correspondence is exact. 

Finally, we have presented experimental evidence for the U{5) — SU{3) transition (spheri- 
cal to axially deformed) in odd-proton nuclei, eiPm., esEu and 65 Tb, and performed a realistic 
calculation which accounts well for the experimental level structure. This is in spite of the 
fact that the structure of these nuclei is rather complex and cannot be simply described 
either by the rotational or by the vibrational model. 

Our investigation of the effects of a single fermion on the phase transitions of a boson 
condensate has been done for the rather complex case of bosons with angular momentum 
L = 0,2 and fermion with angular momentum j = 11/2. It applies equally well to simpler 
cases, for example to the case of a spin j = 1/2 particle immersed in a bath of spinless 
bosons L = 0. As such, the method discussed here can be used in a variety of fields, ranging 
from molecules to atomic condensates, from nuclei to mesoscopic systems. 
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